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Abstract 

Lipid membrane with freely exposed edge is regarded as smooth surface with curved boundary. 
Exterior differential forms are introduced to describe the surface and the boundary curve. The total 
free energy is defined as the sum of Helfrich's free energy and the surface and line tension energy. 
The equilibrium equation and boundary conditions of the membrane are derived by taking the 
variation of the total free energy. These equations can also be applied to the membrane with several 
freely exposed edges. Analytical and numerical solutions to these equations are obtained under the 
axisymmetric condition. The numerical results can be used to explain recent experimental results 
obtained by Saitoh et al. [Proc. Natl. Acad. Sci. 95, 1026 (1998)]. 
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I. INTRODUCTION 



Theoretical study on shapes of closed lipid membranes made great progress two decades 
ago. The shape equation of closed membranes was obtained in 1987 [1], with which the 
biconcave discoidal shape of the red cell was naturally explained Q], and a ratio of \pl of 
he two radii of a torus vesicle membrane was predicted |3[ and confirmed by experiment 



4] 



During the formation process of the cell, either material will be added to the edge or 
the edge will heal itself so as to form closed structure. There are also metastable cup-like 

n 

equilibrium shapes of lipid membranes with free edges p. Recently, opening- up process 
of liposomal membranes by talin [f| Q] has also been observed gives rise to the interest of 
studying the equilibrium equation and boundary conditions of lipid membranes with free 
exposed edges. Capovilla et al. first study this problem and give the equilibrium equation 
and boundary conditions N. They also discuss the mechanical meaning of these equations 

The study of these cup-like structures enables us to understand the assembly process of 
vesicles. Julicher et al. suggest that a line tension can be associated with a domain boundary 
between two different phases of an inhomogeneous vesicle and leads to the budding [1Q ! J . For 
simplicity, however, we will restrict our discussion on open homogenous vesicles. 

In this paper, a lipid membrane with freely exposed edge is regarded as a differentiable 
surface with a boundary curve. Exterior differential forms are introduced to describe the 
surface and the curve. The total free energy is defined as the sum of Helfrich's free energy and 
the surface and line tension energy. The equilibrium equation and the boundary conditions 
of the membrane are derived from the variation of the total free energy. These equations 
can also be applied to the membrane with several freely exposed edges. This is another way 
to obtain the results of Capovilla et al. Some solutions to the equations are obtained and 
the corresponding shapes are shown. They can be used to explain some known experimental 
results 

This paper is organized as follows: In Sec|TIJ we retrospect briefly the surface theory 
expressed by exterior differential forms. In Sec lIII[ we introduce some basic properties of 
Hodge star *. In Sec JIVl we construct the variational theory of the surface and give some 
useful formulas. In SecEJ we derive the equilibrium equation and boundary conditions of 
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the membrane from the variation of the total free energy. In Sec JVI| we suggest some special 
solutions to the equations and show their corresponding shapes. In Sec lVHl we put forward a 
numerical scheme to give some axisymmetric solutions as well as their corresponding shapes 
to explain some experimental results. In Sec JVllll we give a brief conclusion and prospect 
the challenging work. 



II. SURFACE THEORY EXPRESSED BY EXTERIOR DIFFERENTIAL FORMS 

In this section, we retrospect briefly the surface theory expressed by exterior differential 



forms. The details can be found in Ref. 
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We regard a membrane with freely exposed edge as a differentiable and orientational 
surface with a boundary curve C, as shown in Fig^ At every point on the surface, we can 
choose an orthogonal frame e 1; e 2 , e 3 with • ej = 5ij and e 3 being the normal vector. For 
a point in curve C, e± is the tangent vector of C . 

An infinitesimal tangent vector of the surface is defined as 

dr = u x ei + uo 2 e 2 , (1) 

where d is an exterior differential operator, and iO\, u 2 are 1-differential forms. Moreover, we 
define 

de, = uJijej, (2) 

where Uij satisfies = — uiji because of ■ e.,- = dij. 

With dd = and d{u)i Au 2 ) = dui A u 2 — uji A du 2 , we have 

duji = cu 12 A co 2 ; du 2 = uj 2 i t\ uji] uji A uj 13 + cu 2 A cu 23 = 0; (3) 

and 

duij = io ik A u kj = 1, 2, 3), (4) 

where the symbol "A" represents the exterior product on which the most excellent expatia- 
tion may be the Ref. 12 1. 



Eq.© and Cartan lemma imply that: 

u 13 = auji + bcu 2 , cu 23 = bui + ao 2 . (5) 
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Therefore, we have 



area element : dA = u)\ A lu 2 , (6) 
first fundamental form: / = dr ■ dr = uj\ + cu^, (7) 



second fundamental form: 17 = — dr ■ de^ = auj 2 + 2buJiUJ2 + cui 2 , (8) 

„ a + c . . 

mean curvature: H = — - — , (9) 

Gaussian curvature: K = ac — b 2 . (10) 



III. HODGE STAR * 

In this part, we briefly introduce basic properties rather than the exact definition of 
Hodge star * |l3j because we just use these properties in the following contents. 

If g, h are functions defined on 2D smooth surface M, then the following formulas are 
valid: 

*/ = MAu 2 ; (11) 

*df = -/2W1 + fiu 2 , if df = f x u x + f 2 uj 2 ] (12) 

d*df = \7 2 f, V 2 is the Laplace-Beltrami operator. (13) 
We can easily prove that 

// (fd*dg-gd*df)= f {f*dg-g*df) (14) 

J J M JdM 

through Stokes's theorem and integration by parts. 

IV. VARIATIONAL THEORY OF THE SURFACE 

the variation of the surface is defined as: 

5r = Q 2 e 2 + n 3 e 3 , 

where the variation along ei is unnecessary because it gives only an identity. Furthermore, 
let 



Operators d and 5 are independent, thus dS = Sd. d5r = 5dr implies that: 

8uJ\ = ^2^21+^3^31—^2^21, (17) 

5u 2 = dVl 2 + ^3^32 - ^1^12, (18) 

G?f2 3 = fii3Co>i + fi 2 3^2 - ^2^23- (19) 

Furthermore, d5ei = 5ctej implies that: 

5uJij = dVtij + VtikUkj — LUik^kj- (20) 

It is necessary to point out that the properties of the operator 5 are exactly similar to 
those of the ordinary differential. 

V. EQUILIBRIUM EQUATION OF THE MEMBRANE AND BOUNDARY CON- 
DITIONS 

The total free energy F of a membrane with an edge is defined as the sum of Helfrich's 
free energy [l^. \i\ Fh = jj[^f{2H + c ) 2 + kK]dA and the surface and line tension energy 
F s i = A ff dA + 7 § c ds. Here k c , k, c , A and 7 are constants. With the arc-length 
parameter ds — U\, the geodesic curvature k g = uj\ 2 /ds on C and the Gauss-Bonnet formula 
J J KdA = 2tt — § c kgds, the total free energy and its variation are given 

F = JJ[j(2H + c ) 2 + \]w! A u 2 + 7 j> 0J1 - k j> uj 12 + 2nk, (21) 

and 

5F = k c JJ (2H + c^H)^ Au 2 + JJ[j{2H + c ) 2 + A]<f(wi A cu 2 ) 

+ 7 ® 8u>i — k(p 5u>i 2 , (22) 
Jc Jc 

respectively. From Eqs. (fTTjl and (JUJ), we can easily obtain: 

8{uj\ A u) 2 ) = 8uj\ A uj 2 + uj\ A 5uj 2 = —d{VL 2 L)i) — (2H)Q 3 uJi A oj 2 . (23) 

Eqs.©, (HH, (HHI) and flUJ) lead to: 

8{2H)uj\ A uj 2 = 8(a + c)u>i A u 2 

= 2(2H 2 - lT)^i A oo 2 + d(n n uj 2 - fi 23 ^i) 

+ aVL 2 dui — bdQ 2 A u 2 + bQ 2 du 2 + cdQ 2 A u)\. (24) 



Thus we have: 

5F = k c JJ (2H + c ) [2(2H 2 - #)fi 3 ^i Aco 2 + d(n 13 cu 2 - ft 23 ^i) 
+aVL 2 duJi — bdVL 2 /\uj 2 + bQ 2 du 2 + cdQ 2 A uji\ 
+ JJ (^(2H + c ) 2 + A) [-d(O^i) - (2H)Q 3 u 1 A u 2 ] 

+7 (h {tt 2 uj 21 + fi 3 co> 3 i - u 2 {l 2 i] -k(t [dil u + Q 13 uj 32 - u;i 3 f2 32 ]. (25) 
Jc Jc 

If Q 2 = 0, then dQ 3 = + ^23^2, *dQ 3 = —Q 23 uji + ^13^2. On curve C, u 2 = 0, 

^31 = —au>i, u> 32 = —bui, fi 3 |c = ^3c- Thus Eq. (f23J) is reduced to 

5F = [[ [k c {2H + c ) {2H 2 - c H - 2K) - 2A#]ft 3 cu 1 A uo 2 



+ k c JJ (2H + c )d * dtt 3 - 7 J auiVt 3 + k J {bVt V3 - aVt 23 )uoi (26) 
In terms of Eqs. (jl3|) and (fT^j) . we have: 

JJ{2H + c )d * dtt 3 = J (2H + c ) * dtt 3 - J Vt 3 * d{2H + c ) + Jj Vt 3 V 2 {2H + c )^i A co 2 . 

Using integration by parts and Stokes's theorem, we arrive at § c bVLi 3 u)i = § c bdQ 3C = 
— § c VL 3C db. Thus 



SF = jj [k c {2H + c ){2H 2 - c H - 2K) - 2XH + k c V 2 {2H + c )]fi 3 ^i A uj 2 

- (b [k c (2H + c ) + A;a]fi 2 3^i - f ^3c[k c * d(2H + c ) + + kdb]. (27) 
Jc Jc 

It follows that 

k c (2H + c )(2H 2 - c H - 2K) - 2XH + k c V 2 (2H) = 0, (28) 

[k c {2H + co) + ka]\ c = Q, (29) 

[k c *d(2H)+- f auj 1 + kdb}\ c = 0. (30) 

The mechanical meanings of the above three equations are: Eq. (j28j) is the equilibrium 
equation of the membrane; Eq. (j29|) is the moment equilibrium equation of points on C 
around the axis ex" and Eq. (j30|) is the force equilibrium equation of points on C along the 
direction of e 3 jsl 0]. It is not surprising that Eq.(|29p contains the factor k because it is 
related to the bend energy in Helfrich's free energy. However, it is difficult to understand 
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why k is also included in Eq.(j30J). In fact, the term kdb in Eq. (}3T)j) represents the shear stress 
which also contributes to the bend energy in Helfrich's free energy. 

In fact, a = k n and b = r g are the normal curvature and the geodesic torsion of curve C, 
respectively, and *d(2H) = -e 2 • V(2H)u x . Thus Eqs.flU and (EOD become 



[k c {2H + c ) + kk n ]\ c = 0, (31) 

= 0, (32) 



-dr 

[-k c e 2 ■ V(2F) + 7 fc„ + k^ 



c 



respectively. 

If ^3 = 0, then dQ 3 = + ^23^2 — ^2^23 = (^13 — bQ 2 )ui + (^23 — cVL 2 )uj 2 = 0. It 

leads to f2 13 = bQ 2 and f2 2 3 = cQ 2 . 



5F = k c J J (2H + c ) [aVL 2 duji — bdVt 2 A co> 2 + bVt 2 du 2 + cdVL 2 A u)\ + d{Vti 3 uj 2 — 023^1)] 

[^(2H + c ) 2 + X][-d(n 2 iu 1 )]+ 1 ^ n 2 uj 21 -k(£ KQ^. (33) 

Otherwise, u 13 = auji + bu 2 implies that: adu>i + db Au 2 + 2bduj 2 — cdu>i = — daAuJi. Thus 

aVL 2 duji — bdfl 2 A uj 2 + bVL 2 du) 2 + cdQ 2 A ui\ + d(fli 3 u 2 — Q 23 uJi) 
= — d(a + c) A Q 2 ui = —d{2H + Co) A (co is a constant). (34) 



Therefore 



5F = /[-y(2# + co) 2 fi 2 ^i -kKn 2 u x - Afi 2 u;i+7fi 2 W2i]. (35) 
J c 



It follows that: 

= 0, (36) 



[^H + cof + kK + X + ^g 



c 

because of io 2 \ = —k g uJi on C. This equation is the force equilibrium equation of points on 
C along the direction of e 2 0, E| . 

Eqs.pHjl. (J3TJ), (J32*Jl and are the equilibrium equation and boundary conditions of 
the membrane. They correspond to Eqs. (17), (60), (59) and (58) in Ref.jSj, respectively. 
In fact, these equations can be applied to the membrane with several edges also, because in 
above discussion the edge is a general edge. But it is necessary to notice the right direction 
of the edges. We call these equations the basic equations. 
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VI. SPECIAL SOLUTIONS TO BASIC EQUATIONS AND THEIR CORRE- 
SPONDING SHAPES 



In this section, we will give some special solutions to the basic equations together with 
their corresponding shapes. For convenience, we consider the axial symmetric surface with 
axial symmetric edges. Zhou has considered the similar problem in his PhD thesis [16]. 
If expressing the surface in 3-dimensional space as r = {v cosw, v sinw, z(v )} we obtain 



2H=- 



■ sin tp 



fCOS^ 



dv>i 



K 



sin tp cos if) dtp 
v dv 



, V(2tf) 



r2 d l sin tp 



f cos ibf) and V 2 (2H) 



■^f[^os^(^+cos^)], where ib 



arctan 



sec 2 tp dv^ 

^! ], r2 = 9r/dv. Define t as the 



direction of curve C and ri = dr/du. Obviously, t is parallel or antiparallel to ri on curve 
C . Introduce a notation sn, such that sn = +1 if t is parallel to ri, and sn = —1 if not. 



Thus e 2 

k n '- 
to: 



sn 



and e 2 • V{2H) 



sin tp 



T2 

sec tp 



t„ = 0, and k„ = —sn 



-sn< 



cos tp 



cosi/j^(^^ + cosijj-^) on curve C. For curve C, 
. Thus we can reduce Eqs.(J2EJ), ©J, © and © 



. dip , r 1 , sin ib 

— + cosip— - c )[-{ 

v dv 2 v 



,sin if) 



+ cos ib 



dip, 



- A( 



simp 



v 



cosib^-) + k, 
dv 



dv 
cosip d 

v dv 



+ 2 C ° 



.sin^ 



+ cos ib 
v dv 



r . d .smtb .dip-,, 
[v cosip — { h cos^— )\ = U, 



- sin ib 
Co) + k 



' sin^ ,dib 

k c ( h cosip— — 

v dv 



- d smip sm ib 

-snk cos ib— ( ) + 7 

dv v v 



c 



k 2 ,smib. 
2k r v 



o — sin ib cos ib dip 
2 + k --- + A 

v dv 



dv 
= 0, 

c 
0, 

- STry- 



dib. 2 sin if) cos ip dip , 
(37) 



dv 



cosip 



0. 



(38) 
(39) 
(40) 



J c 



In fact, in above four equations only three of them are independent. We usually keep 
Eqs. fjHTj) . (jHS| and (jjQ|) for the axial symmetric surface. For the general case, we conjecture 
that there are also three independent equations among Eqs. (|25j) . (pTT|). (|3^j) and (jHE)- 
Eq.()37p is the same as the equilibrium equation of axisymmetrical closed membranes [l^1l8|. 
In Ref ^j], a large number of numerical solutions to Eq. (j37j) as well as their classifications 
are discussed. 

Next, Let us consider some analytical solutions and their corresponding shapes. We 
merely try to show that these shapes exist, but not to compare with experiments. Therefore, 
we only consider analytical solutions for some specific values of parameters. 
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A. The constant mean curvature surface 



The constant mean curvature surfaces satisfy Eq.(J2EJ) for proper values of k c , cq, K, and 
A. But Eqs.flHD), (E2) and (JUJ) imply 2H + c = 0, k n = 0, and kK + 1 k g = on curve C if 
k c , k, and 7 are nonzero. 

For axial symmetric surfaces, k n = requires sinip = 0. Therefore K = which requires 
kg = 0. Only straight line can satisfy these conditions. It contradicts to the fact that C is 
a closed curve. Therefore, there is no axial symmetric open membrane with constant mean 
curvature. 



B. The central part of a torus 

When A = 0, Co = 0, the condition sin?/> = av + \/2 satisfies Eq. (pIT|) . It corresponds to 
a torus 15]. Eqs. (}38|) and (J40|) determine the position of the edge v e = — , where 

7 (fc c +fc) v /2(2fc2 +4fc c fc+fc 2 ) r- fc 2v/Ti / • + n 4.1, j • • \ 

a = — 1 (2fc +fc)fc • " we ^ kc = k and ^ = (unit: length dimension), 

it leads to 1/a = —1 and v e = 2 ^p- (unit: length dimension). Thus the shape is the 
central part of a torus as shown in FigEl This shape is topologically equivalent to a ring as 
shown in FigC3 



C. A cup 



If we let sinip = according Hu's method [19j Eq.(jS7J) reduce to: 

(v]/ 2 _ i) 1_ \]> ( — y _|_ _i i 1 ( — y 

dv 3 dv 2 dv 2 dv v dv 2 2v dv 

,c 2 2c # A 3^ 2 -2.c^ ,c 2 A 1,* ^ 3 n 
2 u k c 2ir aw 2 k c v z v 2v a 
Now, we will consider the case that ^ = for v — 0. As u — > 0, Eq. (|41|) approaches to 

S~ + «fl — ~^1T W = ®- ^ s s °l u ti° n is ^ = ocx/v + a 2 u + lnu where ai = 0, a 2 and 
a 3 are three constants. If A = and c > 0, we find that ^ = sin-0 = /3(v/v0) + c v ln(t>/t> ) 
satisfies Eq. (|37jl . The shapes of closed membranes corresponding to this solution are fully 
discussed by Liu et flZ.j^nj. Eqs. pSj) and (f4"UJ) determine the position of the edge that satisfies 
tan^(f) = — jk^o ii k = —2k c . If let (3 = 1, v = 1/cq = 1 (unit: length dimension) 
and 7 ^> k c co, we obtain v/v m 1 and its corresponding shape likes a cup as shown in FigEI 
This shape is topologically equivalent to a disk as shown in FigO 



9 



VII. AXISYMMETRICAL NUMERICAL SOLUTIONS 



(2-3 sin 2 ip)if}'v - sin if){\ + cos 2 if}) + [(c 2 + 2X/k c )ip' - {^'f - 2^'"]v 3 



+ [(cq + 2A/fc c )sin'0 - 4c sin ifjip' + 3sin^>(^/) 2 - A cos tfjtfj"}v 2 



k c (c 



sin?/' . -sin?/' 
if) ) — k 



v 

,2 , 



c 



-sin?/' - k , sin ib cosi/j 
kc — k{l + —J h A - swy 



It is extremely difficult to find analytical solutions to Eq.()37p. We attempt to find the 
numerical solutions in this section. But there is a difficulty that sinif>(v) is multi-valued. 
To overcome this obstacle, we use the arc-length as the parameter and express the surface 
as r = {v(s) cosw, v(s) sinw, z(s)}. The geometrical constraint and Eqs.(J2H), (|3*T|) and (|3^|) 
now become: 

v'(s) = cosip(s), z'(s) = sinif>(s), (42) 

(43) 
(44) 

„ ' 2k c ' * ■ - - °- ^ 

We can numerically solve Eqs.([42|) and ()43|l with initial conditions v(0) = 0, ^(O) = 0, 
^'(0) = a and if}"(0) = and then find the edge position through Eqs. (|44|) and (|45j) . The 
shape corresponding to the solution is topologically equivalent to a disk as shown in Fig|2l 
In fact, Eq.(Jinj) can be reduce to a second order differential equation J^, 0,1211, but we still 
use the third order differential equation ()43j) in our numerical scheme. 

In Fig0 we depicts the outline of the cup-like membrane with a wide orifice. The solid line 
corresponds to the numerical result with parameters a = Co = 0.8/im -1 , X/k c = 0.08/im~ 2 , 
"f/k c = 0.20/im" 1 and k/k c = 0.38. The squares come from Fig. Id of Ref.[6j. 

In FiglHl we depicts the outline of the cup-like membrane with a narrow orifice. The 
solid line corresponds to the numerical result with parameters a = cq = 0.86/im _1 , X/k r 
0.26/im -2 , 7//c c = 0.36/im -1 and k/k c = —0.033. The squares come from Fig. 3k of Ref. j 
Obviously, the numerical results agree quite well with the experimental results of Ref. j 



VIII. CONCLUSION 



In above discussion, we introduce exterior differential forms to describe a lipid membrane 
with freely exposed edge. The total free energy is defined as the Helfrich's free energy plus 
the surface and line tension energy. The equilibrium equation and boundary conditions of 
the membrane are derived from the variation of the total free energy. These equations can 
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also be applied to the membrane with several freely exposed edges. A numerical scheme 
to give some axisymmetric solutions and their corresponding shapes do agree with some 
experimental results. 

The method that combines exterior differential forms with the variation of surface is of 
important mathematical significance. It is easy to be generalized to deal with and to simplify 
the difficult variational problems on high-dimensional manifolds. 

Although we give some axisymmetric numerical solutions that agree with experimental 
results obtained by Saitoh et al, up to now, we still cannot find any unsymmetrical solution. 
A large number of unsymmetrical shapes are found in experiments, which will be a challenge 
to the theoretical study. 
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FIG. 1: The surface with an edge C. At every point of the surface, we can construct an orthogonal 
frame ei,e2,e3, where e3 is the normal vector of the surface. For a point on curve C, e\ is the 
tangent vector of C. 




FIG. 2: The central part of a torus. 
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FIG. 3: A ring (left) and a disk (right). 




FIG. 5: The outline of the cup-like membrane with a wide orifice. The solid line is the numerical 
result with parameters a = cq = 0.8/um -1 , X/k c = 0.08/xm -2 , j/k c = 0.20/im -1 and k/k c = 0.38. 

Q 

The squares come from Fig. Id of Ref. 6]. z-axis is the revolving axis and v is the revolving radius. 
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FIG. 6: The outline of the cup- like membrane with a narrow orifice. The solid line is the numerical 
result with parameters a = cq = O.SQfim^ 1 , X/k c = 0.26/um -2 , ^/k c = O^Qfim" 1 and k/k c = 
—0.033. The squares come from Fig. 3k of Ref. j3|. z-axis is the revolving axis and v is the revolving 
radius. 
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